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ABSTRACT 

Idiopathic osteoporosis (lOP) in premenopausal women is a poorly understood entity in which otherwise healthy women have low- 
trauma fracture or very low bone mineral density (BMD). In this study, we applied individual trabeculae segmentation (ITS)-based 
morphological analysis to high-resolution peripheral quantitative computed tomography (HR-pQCT) images of the distal radius and 
distal tibia to gain greater insight into skeletal microarchitecture in premenopausal women with lOP. HR-pQCT scans were performed for 
26 normal control individuals and 31 women with lOP. A cubic subvolume was extracted from the trabecular bone compartment and 
subjected to ITS-based analysis. Three Young's moduli and three shear moduli were calculated by micro-finite element {|xFE) analysis. 
ITS-based morphological analysis of HR-pQCT images detected significantly decreased trabecular plate and rod bone volume fraction 
and number, decreased axial bone volume fraction in the longitudinal axis, increased rod length, and decreased rod-to-rod, plate-to-rod, 
and plate-to-plate junction densities at the distal radius and distal tibia in women with lOP. However, trabecular plate and rod thickness 
did not differ. A more rod-like trabecular microstructure was found in the distal radius, but not in the distal tibia. Most ITS measurements 
contributed significantly to the elastic moduli of trabecular bone independent of bone volume fraction (BV/TV). At a fixed BV/TV, plate- 
like trabeculae contributed positively to the mechanical properties of trabecular bone. The results suggest that ITS-based morphological 
analysis of HR-pQCT images is a sensitive and promising clinical tool for the investigation of trabecular bone microstructure in human 
studies of osteoporosis. © 2010 American Society for Bone and Mineral Research. 
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Introduction 



Osteoporosis, characterized by low bone mass and micro- 
architectural deterioration, is diagnosed most commonly in 
postmenopausal women and elderly men.'^* However, young or 
middle-aged men and premenopausal women can also suffer 
from osteoporosis. Most young individuals with osteoporosis 
have an underlying secondary cause, such as an endocrine or 
metabolic disorder (eg, hypogonadism) or medication exposure 
(eg, glucocorticoids).'^"^ ^' However, 40% to 60% of young men 
and women with fractures and/or low areal bone mineral density 



(aBMD) by dual-energy X-ray absorptiometry (DXA) have no 
discernible secondary cause.'^'^"'^^"^^' Such individuals are 
operationally considered to have idiopathic osteoporosis (lOP). 

Measurement of aBMD by DXA is the current gold standard for 
diagnosing osteoporosis in postmenopausal women.*^^' How- 
ever, the clinical significance of low aBMD measurements in 
premenopausal women is uncertain, particularly in those who do 
not have a history of fragility fractures. Low aBMD measurements 
may be related to small bone size.'^'^'^^' They may also be due to 
low peak bone mass, either genetically determined or related to 
suboptimal bone accrual during adolescence.'^^' Under these 
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circumstances, bone microarchitecture and strength may be 
normal and short-term fracture risk low. Conversely, low aBMD 
may be due to premature bone loss and associated with 
compromised bone microarchitecture and reduced bone 
strength. Newer technologies, capable of imaging the skeleton 
three-dimensionally with higher resolution than DXA, are 
required to determine whether there are specific microarchi- 
tectural abnormalities In premenopausal women with lOP. 

Recently, an In vivo clinical imaging modality, high-resolution 
peripheral quantitative computed tomography (HR-pQCT), has 
been developed to assess bone microstructure (XtremeCT, Sca- 
nco IVledical AG, Bassersdorf, Switzerland). The 3D data sets 
provided by HR-pQCT permit the separate analysis of cancellous 
and cortical bone at peripheral sites (distal tibia and radius).'^" 
The resolution of HR-pQCT allows for the visualization and 
quantitative evaluation of the internal structure of trabecular 
bone. Several clinical studies have demonstrated the ability of 
HR-pQCT to detect age- or disease-related changes in bone 
microarchitecture and to provide additional fracture risk determi- 
nants.'''"^'" In addition, we recently applied HR-pQCT technology 
to a group of premenopausal women with IQP and detected 
abnormal trabecular microstructure and reduced mechanical 
competence.'^^' 

Despite its great potential, the current standard morphological 
analysis of HR-pQCT has several limitations. First, the resulting 
morphological parameters are based on a statistical average of 
the imaged trabecular network rather than on measurements of 
individual trabeculae. Second, the current analysis technique 
evaluates the morphology of trabecular bone at a global level 
without separate analyses for trabecular plates and rods, even 
though there are fundamental differences between these two 
types of structures. It is becoming increasingly apparent that the 
microstructural type of trabeculae (the two types of structures 
are plate-like and rod-like) is critical in determining the strength 
of trabecular bone'^''"^^' because a dramatic change in 
trabeculae from plate-like to rod-like occurs with aging and 
osteoporosis.''^*"^'* Furthermore, structure model index (SMI), a 
parameter used in micro-computed tomography (jaCT) analysis 
to characterize the plate-likeness of the trabecular network,'^"''^' 
is not included in the standard analysis output of HR-pQCT scans. 
Most importantly, several HR-pQCT standard parameters are 
derived rather than measured directly.'^^' For example, the 
trabecular number (Tb.N*) parameter is defined as the inverse of 
the mean distance between the midline of trabeculae. Trabecular 
thickness (Tb.Th) and trabecular spacing (Tb.Sp) are derived 
based on Tb.N* and the derived bone volume fraction (BV/TV^) 
[ie, Tb.Th = BV/TV'^'/Tb.N* and Tb.Sp=(1 - BV/TV^j/Tb.N*] by 
analogy to standard histomorphometry.'^^' Therefore, since 
these parameters are all highly dependent on BV/TV^, they may 
only provide limited additional information. 

In our previous work we studied the relative importance of 
trabecular types (plates and rods) in the architecture and 
mechanical properties of trabecular bone and developed image- 
processing techniques for volumetrically segmenting the 3D 
trabecular bone microstructure as a collection of trabecular 
plates and rods.'^^'^^' We developed individual trabeculae 
segmentation (ITS)-based morphological analysis to assess 
trabecular microstructure based on measurements of each 



individual plate and rod.' ' One clinical application of ITS 

analysis was based on in vivo micro-magnetic resonance (jjlMR) 
images of tibial trabecular bone: ITS-based morphological 
analysis detected significant trabecular microarchitectural 
changes in plates and rods in response to testosterone treatment 
in hypogonadal men.'^^' With the recent success of HR-pQCT in 
assessing microstructure and mechanical competence in 
osteoporosis, there is great potential in the application of ITS 
technology to clinical HR-pQCT images for in vivo and model- 
independent assessment of trabecular microarchitecture with a 
detailed quantification of trabecular types. 

In this study, we applied ITS-based morphological analysis to 
HR-pQCT images of the distal radius and tibia to gain greater 
insight into trabecular bone microstructure changes in otherwise 
healthy premenopausal women with lOP. Then, the contribu- 
tions of ITS-based microstructural parameters of HR-pQCT to 
mechanical properties estimated by micro-finite element (ijlFE) 
analysis of trabecular bone were evaluated. We hypothesized 
that trabecular microarchitecture in lOP would be characterized 
by conversion of trabecular plates into rods and decreased 
trabecular network connectivity and trabecular number, rather 
than uniform trabecular thinning. 

Materials and Methods 

Patient population 

Premenopausal women (aged 20 to 48 years, 26 control and 
31 lOP) were recruited at Columbia University and Creighton 
University. Subjects with lOP were included on the basis of a 
documented history of low-trauma fractures after age 18 and/or 
very low bone mineral density (BMD) measurement (T-score < 
—2.5 or Z-score < —2.0) at the lumbar spine, total hip, or femoral 
neck. Low-trauma fractures were defined as equivalent to a fall 
from a standing height or less, excluding skull or digit fractures. 
Fractures were ascertained by review of radiographs or radio- 
graph reports. To qualify as normal controls, women were 
required to have normal aBMD by DXA (T-score > —1.0) and no 
history of low-trauma fractures. Subjects were excluded if they 
had given birth or had lactated within the past 12 months. 
Patients and controls had a detailed history, physical examina- 
tion, and biochemical evaluation to exclude secondary causes of 
osteoporosis, including disorders causing premenopausal estro- 
gen deficiency, endocrinopathies (eg, hyperthyroidism, Cushing 
syndrome, prolactinoma), anorexia nervosa, bulemia, celiac or 
other gastrointestinal diseases, abnormal mineral metabolism 
(eg, osteomalacia, hyperparathyroidism), vitamin D deficiency, 
hypercalciuria (>300 mg/g of creatinine), and drug exposures (eg, 
glucocorticoids, anticonvulsants, anticoagulants, methotrexate). 

HR-pQCT innages of the distal radius and distal tibia 

HR-pQCT (XtremeCT, Scanco Medical AG) was performed at 
Columbia University for all participants, including those recruited 
at Creighton University. The nondominant forearm and distal 
tibia were immobilized in a carbon-fiber shell and scanned as 
described previously.'^^'^"'^'^* The region of interest was defined 
on a scout view by manual placement of a reference line at the 
endplate of the radius or tibia — the first slice being 9.5 and 
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22.5 mm proximal to the reference line at the radius and tibia, 
respectively. The HR-pQCT measurement included 110 slices 
corresponding to a 9.02-mm section along the axial direction 
with a nominal voxel size of 82 jjim. The European forearm 
phantom was scanned whenever subjects were scanned for 
quality control. 

According to the standard patient evaluation protocol, the 
periosteal surfaces of the radius and tibia were contoured first 
semi-automatically and then by an automated threshold-based 
algorithm to separate the cortical and trabecular compart- 
ments.'^^' The mineralized phase then was thresholded auto- 
matically using a Laplace-Hamming filter followed by a global 
threshold using a fixed value of 40% of the maximal grayscale 
value of the images.'^^' These procedures were performed using 
Scanco software installed on an HP AlphaStation operating in a 
VMS environment (Hewlett-Packard, Palo Alto, CA, USA). After 
each image was thresholded, a 5.7 x 5.7 x 5.7 mm^ cubic subvol- 
ume corresponding to 70 x 70 x 70 voxels was extracted from 
the center of trabecular bone compartment of each image of the 
distal radius; similarly, a 9.0 x 9.0 x 9.0 mm^ subvolume corre- 
sponding tollOxllOxllO voxels was extracted from each 
distal tibia (Fig. 1). Owing to the size difference of the distal radius 
and distal tibia segment, the cubic subvolume of the distal radius 
was smaller than that of the distal tibia. 

Individual trabeculae segmentation (ITS)-based 
morphological analysis of HR-pQCT images 

All the trabecular bone subvolumes of the distal radius and distal 
tibia from the control and lOP subjects were subjected to ITS- 
based morphological analysis. First, a complete volumetric 
decomposition technique was applied to segment the trabecular 
network into individual plates and rods.'^^' Briefly, digital 
topologic analysis (DTA)-based skeletonization*^^' was applied 
first to transform a trabecular bone image into a representation 
composed of surfaces and curves skeleton while preserving the 
topology (ie, connectivity, tunnels, and cavities),'^^'^^' as well as 
the rod and plate morphology of the trabecular microarchi- 
tecture. Then digital topologic classification was applied in which 
each skeletal voxel was uniquely classified as either a surface or a 
curve type.'^^' Using a newly developed iterative reconstruction 
method, each voxel of the original image was classified as 
belonging to either an individual plate or rod (Fig. 1). Based on 
the evaluations of dimension and orientation of each individual 



trabecular plate and rod, as well as junctions of surface and curve 
skeletons, a set of ITS-based morphological parameters was 
derived to quantify plate and rod bone volume fraction (pBV/TV 
and rBV/TV), axial bone volume fraction along the longitudinal 
axis (aBV/TV), plate and rod tissue fraction (pBV/BV and rBV/BV), 
that is, the volume of plate/rod bone tissue divided by the total 
volume of bone tissue, plate and rod number density (pTb.N and 
rTb.N, 1/mm), plate and rod thickness (pTb.Th and rTb.Th, mm), 
plate surface area (pTb.S, mm^), rod length [rJhl, mm), and rod- 
rod, plate-rod, and plate-plate junction density (R-R, P-R, and P-P 
Junc.D, l/mm^).The definition of these ITS measurements can be 
found in the Glossary at the end of this article. Detailed methods 
of the complete volumetric decomposition technique and ITS- 
based measurements can be found in our recent publication.'^^' 

Finite element analysis for HR-pQCT images 

Each subvolume of HR-pQCT image of the distal radius and distal 
tibia was converted to a |jlFE model by directly converting bone 
voxels to 8-node elastic brick elements with an element size of 
82 X 82 X 82 (jum^. Using a customized element-by-element 
preconditioned conjugate gradient solver,''"" six |jlFE analyses 
were performed for each model, representing three uniaxial 
compression tests along three imaging axes (x, y, and z) and 
three uniaxial shear tests.''*^' The general anisotropic stiffness 
matrix was first determined based on the results from the 
preceding analyses. A new coordinate system of orthotropic axes 
(Xi, X2, and X3) representing the best orthotropic symmetry was 
then calculated using a numerical optimization algorithm, 
Powell's method,'*^' to minimize an orthotropy objective 
function.''"' Transformation of the anisotropic matrix to a new 
coordinate system yielded the full orthotropic stiffness tensor.'''^' 
The elastic constants and stiffness matrix were sorted such that 
£11 was in the direction of the lowest axial modulus and £33 was 
in the direction of the highest axial modulus (longitudinal axis). 
The elastic moduli (three Young's moduli, £ii<£22<£33, and 
three shear moduli, G23, 631, and 612) were then derived from the 
orthotropic stiffness tensor. All the jjiFE analyses for subvolumes 
were implemented on a Dell XPS PC workstation (Dell, Inc., 
Round Rock, TX, USA). 

Statistical analysis 

A nonparametric statistical test, a two-tailed Mann-Whitney U 
test, was performed to test for significant differences in the ITS 
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measurements and the six elastic moduli between the control 
and lOP groups, as well as between the low BMD and fracture 
subjects in the lOP group. To investigate the contributions of the 
ITS measurements to the mechanical properties of trabecular 
bone, data from both the control and lOP groups were pooled. 
Each of the ITS measurements was correlated individually to the 
elastic moduli of trabecular bone at both the distal radius and 
distal tibia. Next, to differentiate the contributions of the ITS 
measurements to the mechanical properties of trabecular bone 
from that of BV/TV, partial correlation analyses between elastic 
moduli and ITS measurements with the effect of BV/TV removed 
were performed. A forward stepwise multiple linear regression 
was also performed to predict the elastic moduli by the ITS 
measurements. From the ITS-based parameters, pBVATV, rBV/TV, 
aBV/TV, pBV/BV, pTb.N, rTb.N, pTb.Th, rTb.Th, pTb.S, rTbl, P-P 
Junc.D, P-R Junc.D, and R-R Junc.D were used as the independent 
variables. Finally, correlations between each of the ITS 
measurements were calculated for both the distal radius and 
tibia. All the statistical analyses were performed using 
KaleidaGraph 3.6 software (Synergy Software, Reading, PA, 
USA) or SPSS 13.0 software (SPSS, Inc., Chicago, IL, USA). 

Results 



The study population consisted of 26 control individuals and 31 
lOP subjects Twenty-one of these subjects were included 
because they experienced low-trauma adult fractures whether 
or not they had low aBMD (fracture group), and the other 1 0 had 
low aBMD but no adult low-trauma fractures (low BMD group). 
The mean ages of the control and lOP subjects were 34.3 ± 8.7 
and 37.8 ± 7.8 years, respectively. The age, height, weight, and 
body mass index (BMI) did not differ between groups (p > .05). 
The mean Z-scores at the lumbar spine, total hip, femoral neck. 



and distal radius were 0.61 ±0.82, 0.28 ±0.63, 0.1 2 ±0.75, and 
0.71 ±0.80 in the control group and -1.42± 1.16, -1.00± 1.30, 
-1.35 ± 1.22, and 0.1 8 ±0.82 in the lOP group. Although most 
lOP subjects were included on the basis of low-trauma fractures 
regardless of aBMD, aBMD was lower at all sites in lOP subjects 
than in the control individuals (p < .05). Detailed characteristics 
of the study population were reported previously.'^^' 

The mean and standard deviation of the HR-pQCT image- 
based ITS and mechanical measurements of the control and lOP 
subjects are presented in Table 1 and Fig. 2. At the distal radius, 
plate bone volume fraction (pBVATV), rod bone volume fraction 
(rBV/TV), and axial bone volume fraction (aBV/TV) were 42%, 
23%, and 30% lower, and plate number (pTb.N) and rod number 
(rTb.N) were 20% and 10% lower in women with lOP. In addition, 
plate tissue fraction (pBV/BV) was 27% lower in women with lOP. 
In contrast, rod length (rTbi) was 1 0% higher and plate surface 
area (pTb.S) did not differ in lOP subjects. Rod-to-rod, plate-to- 
rod, and plate-to-plate junction densities (R-R, P-R, and P-P 
Junc.D) were 26%, 42%, and 45% lower in women with lOP. 
Similar results were found at the distal tibia. pBV/TV was 23% 
lower but not statistically significant (p = .l), and rBV/TV and 
aBV/TV were 31 % and 1 8% lower in women with lOP. pTb.N and 
rTb.N were both 1 6% lower in women with lOP. However, pBV/BV 
did not differ. pTb.S and rJhl were 10% and 7% higher, and R-R, 
R-P, and P-P Junc.D values were 37%, 37%, and 34% lower in 
women with lOP. Neither pTb.Th nor rTb.Th differed in women 
with lOP at the distal radius or distal tibia (Table 1). At both the 
distal radius and the tibia, the six elastic moduli of the lOP group 
were significantly lower than those of the control group (Fig. 2). 
No statistical difference was detected in BV/TV, ITS measure- 
ments, or elastic moduli between fracture and low BMD subjects 
in the lOP group. 

At the distal radius, the ITS measurements, including pBV/TV, 
rBV/TV, aBV/TV, pBV/BV, pTb.N, rTb.N, pTb.Th, rTbJh, and R-R, P- 



Table 1. ITSMeasurements(Mean±SD)ofHR-pQCT Trabecular Bone Images at the Distal Radius and the Distal Tibia of Control and lOP 
Groups 







Distal radius 






Distal tibia 




Control (n — 26) 


lOP (fi-31) 


P 


Control (n — 26) 


lOP (fi-31) 


P 


pBV/TV 


0.053 ± 0.026 


0.031 ±0.028 


.0013 


0.043 ±0.01 9 


0.033 ± 0.027 


.0232 


rBV/TV 


0.171 ±0.040 


0.1 32 ±0.039 


.0002 


0.147 ±0.034 


0.101 ±0.037 


.0001 


aBV/TV 


0.077 ±0.025 


0.054 ±0.029 


.0008 


0.067 ±0.01 5 


0.055 ±0.024 


.0193 


pBV/BV 


0.22 ± 0.09 


0.1 6 ±0.09 


.006 


0.22 ± 0.08 


0.23 ±0.11 


.94 


rBV/BV 


0.78 ±0.09 


0.84 ±0.09 


.006 


0.78 ±0.08 


0.77 ±0.11 


.94 


pTb.N (1/mm) 


1.21 ±0.25 


0.97 ± 0.27 


.0007 


1.13±0.15 


0.96 ±0.25 


.0023 


rTb.N (1/mm) 


1 .86 ± 0.22 


1.67 ±0.21 


.0002 


1.76 ±0.1 9 


1 .48 ± 0.27 


.0001 


pTb.Th (mm) 


0.1 92 ±0.007 


0.1 88 ±0.011 


.20 


0.200 ±0.008 


0.202 ±0.014 


.82 


rTb.Th (mm) 


0.21 3 ±0.007 


0.21 3 ±0.008 


.76 


0.21 5 ±0.008 


0.218 ±0.008 


.22 


pTb.S (mm'^) 


0.138±0.107 


0.1 42 ±0.1 52 


.14 


0.1 38 ±0.1 30 


0.1 53 ±0.289 


.02 


rThl (mm) 


0.674 ± 0.048 


0.734 ±0.066 


.0002 


0.689 ±0.038 


0.736 ±0.046 


.0002 


R-R JuncD (1/mm^) 


3.23 ±0.90 


2.39 ±0.88 


.0008 


2.66 ±0.93 


1.67 ±0.78 


.0001 


P-R Junc.D (1/mm^) 


2.85 ±1.10 


1.64± 1.16 


.0004 


2.23 ± 0.74 


1 .40 ± 0.83 


.0003 


P-P Junc.D (1/mm^) 


1 .23 ± 0.54 


0.67 ±0.55 


.0005 


0.96 ± 0.36 


0.63 ± 0.45 


.0014 



Note: p values were calculated by two-tailed Mann-Whitney U tests for comparison between groups. SD = standard deviation. Please refer to the Glossary 
for definitions of each ITS-based parameter. 
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Fig. 2. Comparisons of estimated elastic moduli between control and lOP groups at (A) the distal radius and (B) the distal tibia. The asterisk indicates a 
significant difference between control and lOP groups (two-tailed Mann-WhItney U test, p < .05). 



R, and P-P Junc.D, were significantly and positively correlated 
with BV/TV and the six elastic moduli of trabecular bone, whereas 
rBV/BV and rXbl were negatively correlated with BV/TV and the 
elastic moduli (p<.001; Table 2). pTb.S was weakly and neg- 
atively correlated with BV/TV but not statistically correlated with 
any elastic modulus (p > .05). The results from the distal tibia 
were similar to those from the distal radius with a few exceptions: 
rTb.Th had no contribution to any elastic modulus, and neither 
pTb.Th nor rTb.Th were correlated with BV/TV (Table 2). For both 
the distal radius and the distal tibia, the strongest correlations 
were between the axial modulus £33 and pBV/TV, aBV/TV, pTb.N, 
P-R Junc.D, and P-P Junc.D. 

Partial correlation analysis with the effect of BV/TV removed 
revealed that each ITS measurement contributed independently 
to the elastic moduli of trabecular bone. At the distal radius, pBV/ 
TV, aBV/TV, pBV/BV, pTb.N, rTb.Th, rTbl, P-R Junc.D, and P-P 
Junc.D had significant and positive partial correlations with most 
elastic moduli of trabecular bone. rBV/TV, rBV/BV, rTb.N, and R-R 
Junc.D showed significant and negative partial correlations with 
most elastic moduli of trabecular bone. pTb.Th had positive and 
significant partial correlations with £33 but no correlation with 
the other moduli. At the distal tibia, pBV/TV, aBV/TV, pBV/BV, 
pTb.N, pTb.Th, P-R Junc.D, and P-P Junc.D had significant and 
positive partial correlations with most elastic moduli, whereas 
rBV/TV, rBV/BV, rTb.N, and R-R Junc.D had negative correlations. 
rTb.^ had significant and positive partial correlation with £33 but 
not with the other moduli (Table 2). Furthermore, multiple linear 
regression analysis revealed that each ITS-based microstructural 
parameter independently contributed to one or several elastic 
moduli of trabecular bone (Table 3). 

The correlations between ITS-based parameters are listed in 
Table 4. In general, plate-related parameters, including pBV/TV, 
pBV/BV, pTb.N, pTb.Th, P-R Junc.D, and P-P Junc.D, as well as aBV/ 
TV, were significantly and positively correlated with each other at 
both the distal radius and the tibia. These parameters also 
significantly and negatively correlated with rBV/BV and rTb.^. 

Discussion 



In this study, ITS-based morphologic analysis, which can provide 
model-independent and trabecular type-specific evaluations 



of trabecular bone microarchitecture, was applied to in vivo 
HR-pQCT images of the distal radius and tibia from premeno- 
pausal women with and without lOP. We found that pre- 
menopausal women with lOP had fewer longitudinally aligned 
trabeculae, fewer trabecular plates, fewer and longer trabecular 
rods, and decreased connectivity between rods and plates at 
both the distal radius and tibia. While both plate and rod bone 
volume fractions were significantly lower, trabecular plate and 
rod thickness did not differ in women with lOP at either the distal 
radius or distal tibia. These results suggest that there has been a 
loss of trabecular plates and rods in women with lOP, resulting in 
a less connected, less longitudinally aligned, and more widely 
separated trabecular network and reduced bone mass. Further- 
more, the plate tissue fraction at the distal radius was 
significantly lower in women with lOP compared with controls, 
indicating trabecular conversion from plates to rods in the 
trabecular network. However, this change in plate tissue fraction 
was not found at the distal tibia. Within the lOP group, there was 
no statistical difference between fracture and low BIVID subjects 
in any of the ITS-based parameters or elastic moduli. It is 
noteworthy to mention that the trabecular bone microarchi- 
tecture and mechanical properties in women with low BMD who 
had never had an adult low-trauma fracture were as severely 
affected as those who had had low-trauma fractures. 

Previous retrospective transiliac biopsy studies of young 
adults with unexplained fractures have yielded conflicting data 
on bone microarchitecture.'""^^' In a recent report,''^^' we exa- 
mined in vivo microstructure measurements using the standard 
patient analysis protocol of HR-pQCT in the same populations of 
control and lOP subjects as the current study. Severe trabecular 
microarchitecture deterioration of the distal radius and distal 
tibia was detected in premenopausal women with lOP — redu- 
ced trabecular number and thickness and increased trabecular 
spacing and network inhomogeneity. The data reported in this 
study confirm the major findings we recently reported based on 
the standard patient analysis but also refine and extend them.''^^' 

High-resolution 3D images by HR-pQCT can reflect most 
trabecular microarchitecture when compared with gold stand- 
ard |jiCT images'"''*'; however, their potential is currently limited 
by the built-in microstructural analysis technique. Measurements 
such as Tb.Th and Tb.Sp are model-dependent. Furthermore, no 
measurements of trabecular type and trabecular network con- 
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Table 3. Correlation Coefficient and Independent Predictors of Multilinear Regression for the Predictions of Elastic Moduli of Trabecular 
Bone by the ITS Measurements 







Distal radius 




Distal tibia 


Mechanical 


Adjusted 




Adjusted 




properties 


r 


Independent predictors 


r 


Independent predictors 


c 

Eu 


0.96 


P-R Junc.D, rTD.Tn, pTb.N, rTb.*, pTb.Tn 


0.93 


K-K Junc.D, riD. m, K-K Junc.U, pdv/dv 


E22 


0.98 


rJb.Th, pBV/BV, pBV/Tv, rTb.£, rlb.N 


0.95 


P-P Junc.D, rBV/TV, rTb.N, P-R Junc.D 


E33 


0.99 


P-P Junc.D, rTb.Th, aBv/TV, pTb.N, rTb.t; 


0.99 


P-P Junc.D, aBV/TV, rTb.Th, R-R Junc.D, 










pTb.Th,* pTb.N,** pTb.S*** 


G23 


0.99 


P-P Junc.D,* rTb.Th,* pBV/BV,**** 


0.96 


P-P JuncD* 






pBV/TV,*** rTb.£,* pTb.N* 






G31 


0.98 


P-P Junc.D,* rBV/TV,* rlb.N,* rTb.£** 


0.96 


rTbJh,* R-R Junc.D,* pBV/TV,* rTb.N** 


G12 


0.98 


P-R Junc.D,*** rTb.Th,*** pTb.N,* 


0.94 


P-R JuncD* 






rVoi;* pBV/TV,* rBV/TV,* rTb.N** 







Note. The order of predictors Is given in the order that predictors were entered in the forward stepwise regression model. The significance of each 
independent predictor Is Indicated by Its p value (*p < .001 ; **p < .01 ; ***p < .05, ****p < .1 ). Please refer to the Glossary for definitions of each ITS-based 
parameter. 



nectivity are provided by HR-pQCT standard analysis. In this 

study, we used a model-Independent ITS technique to examine 
the differences In trabecular microstructure at the individual 
trabecula level. The ITS technique confirmed our previous 
findings using routine analysis of HR-pQCT scans, specifically that 
lower trabecular number and a more separated trabecular 
network underlie lower bone mass in premenopausal women 
with lOP. However, ITS-based analysis also provided additional 
detail about the trabecular structure In lOP, revealing a lower 
degree of trabecular alignment and loss of trabecular con- 
nectivity between trabecular rods and rods, rods and plates, and 



plates and plates in women with lOP. In contrast to standard HR- 
pQCT analysis measurements, which showed a significantly 
reduced Tb.Th, ITS analysis indicates that plate and rod thickness 
remained unchanged. There are two possible explanations for 
the contradiction. First, by HR-pQCT standard analysis, Tb.Th Is a 
derived parameter, based on Tb.N*, calculated as the Inverse of 
the distance between trabeculae midlines, and BV/TV^, a derived 
bone volume fraction assuming a density of fully mineralized 
bone tissue of 1200mg HA/cm^. Therefore, the accuracy of the 
Tb.Th measurement by HR-pQCT standard analysis is subject to 
several assumptions. The variation may be even greater in 



Table 4. Correlation Coefficient of Linear Regression Between Each of the ITS Measurements of the HR-pQCT Images of Trabecular Bone 
at Both the Distal Radius and the Distal Tibia 




pBV/ rBV/ aBV/ pBV/ 
TV TV TV BV 

pBV/TV ■■^^^^■0.96* 

rBV/TV 0.40** ^^Ho.45* 

aBV/TV 0.88* 0.57* 

pBV/BV 0.83* NS 0.61* 

rBV/BV -0.83* NS -0.61* -1.00* 

pTb.N 0.90* 0.67* 0.91* 0.61* 

rTb.N 0.39** 0.97* 0.56* NS 

pTb.Th 0.63* NS 0.50* 0.82* 

rTb.Th NS -0.31*** NS 0.32*** 

pTb.S NS -0.61* NS 0.57* 

llbl -0.53* -0.85* -0.54* NS 

R-R Junc.D NS 0.96* 0.36** NS 

P-R Junc.D 0.84* 0.78* 0.83* 0.48* 

P-P Junc.D 0.93* 0.66* 0.87* 0.62* 



rBV/ 
BV 

-0.96* 
-0.45* 



pTb.N rTb.N pTb.Th 

Distal radius 

0.94* 0.55* 0.71* 

0.79* 0.98* 0.50* 

).97* 0.67* 0.71* 

0.41** 0.71* 

-0.41** -0.71* 



-0.61* 

NS 0.70* 

-0.82* 0.44* NS 

-0.32*** NS -0.42** 0.42** 

-0.57* NS -0.71* 0.59* 

NS -0.71* -0.83* NS 

NS 0.50* 0.93* NS 

-0.48* 0.94* 0.75* 0.30*** 

-0.62* 0.95* 0.64* 0.41** 
Distal tibia 



rTb.Th pTb.S 




0.39** 

0.42** 0.34*** 

-0.51* -0.63* 

NS NS 

NS NS 



R-R P-R P-P 
rTb.^ JuncD Junc.D JuncD 



0.81* 0.73* 

0.95* o.9r- 

0.83* 0.8* 

-0.83* -0.89* 

0.97* 0.97* 

0.77* 0.6H 



NS 


-0.72* 


0.45* 


0.47* 


-0.92* 


0.97* 


NS 


-0.81* 


0.57* 


NS 


-0.62* 


0.28*** 


NS 


0.62* 


-0.28*** 


NS 


-0.88* 


0.67* 


0.46* 


-0.93* 


0.97* 


NS 


-0.54* 


0.34** 






_ NS 



0.38* 
NS 



52* 

-0.85* 

-0.82* 0.66* 
-0.74* 0.51* 0.98* 



0.45*1 
N^ 



-0.87* -0.82* 



Wote: Data were combined from both control and lOP groups. Data for the distal radius are highlighted In light gray. 
*p < .001 ; **p < .01 ; ***p < .05. NS: p > .05. Please refer to the Glossary for definitions of each ITS-based parameter. 
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pathologic conditions in which tissue density is altered. Second, 
however, compared with mean trabecular thickness (100 to 
300 (Am), HR-pQCT image voxel size (82 [xm) may not be 
sufficient for ITS analysis to detect differences in trabecular 
thickness, especially if they are subtle. 

The three Young's moduli and three shear moduli were all 
significantly lower in the lOP group than in the control group. 
Thus mechanical properties of trabecular bone in premenopau- 
sal women with lOP were severely compromised. IVlost ITS 
measurements were good indicators of the elastic moduli of 
trabecular bone. In general, linear correlations between the 
elastic moduli and ITS measurements suggested that pBV/TV, 
rBV/TV, aBV/TV, pBV/BV, pTb.N, rTb.N, pTb.Th, R-R Junc.D, P-R 
Junc.D, and P-P Junc.D correlate positively, whereas rBV/BV and 
rJb.£ correlate negatively with the mechanical properties of 
trabecular bone. Although most ITS measurements were also 
significantly associated with BV/TV (Table 2), partial correlation 
results showed that at a fixed BV/TV, plate-associated parameters 
such as pBV/TV, aBV/TV, pBV/BV, pTb.N, pTb.Th, P-R Junc.D, and 
P-P Junc.D contribute positively to the mechanical properties of 
trabecular bone at both the distal radius and the distal tibia. In 
contrast, rod-associated parameters such as rBV/TV, rBV/BV, 
rTb.N, and R-R Junc.D correlate negatively with trabecular bone's 
mechanical properties. This is consistent with our findings based 
on |xCT images of human trabecular bone that trabecular 
plates play an important role in determining the mechanical 
properties of bone and that a conversion of trabecular network 
from plates to rods impairs bone strength.*^^'^*' 

There were also interdependencies between most of the ITS 
measurements (Table 4). For example, an increase of pBV/TV is 
associated with an increase of pBV/BV, pTb.N, P-P Junc.D, pTb.Th, 
and aBV/TV, as well as a decrease of rBV/BV and rTb.^. However, 
although associated with each other, each ITS measurement 
independently contributes to one or more elastic moduli, as 
demonstrated by multiple linear regressions of the elastic moduli 
with ITS-based parameters (Table 3). The combination of ITS 
measurements showed excellent predictions of the mechanical 
properties of trabecular bone, with correlation coefficient r 
ranging from 0.96 to 0.99 at the distal radius and 0.93 to 0.99 at 
the distal tibia. 

In this study, the elastic moduli of trabecular bone were 
estimated by HR-pQCT image-based |jlFE analysis based on a 
subvolume of trabecular region at the distal radius and the distal 
tibia. This (xFE analysis technique has been thoroughly validated 
on HR-pQCT images with the registered |xCT images as a gold 
standard.'*^' Correlations between HR-pQCT- and jjiCT-based 
elastic moduli of trabecular subvolumes were highly significant, 
with > 0.90.''*^' The accuracy of ITS measurements was also 
investigated in a preliminary validation study.'"*^' Most ITS 
measurements of HR-pQCT images correlated significantly with 
those of the registered [xCT images, and plate-related 
parameters yielded better correlations than rod-related para- 
meters. This suggests that the partial-volume effect in the low- 
resolution regimen has a greater influence on trabecular rods 
than plates. However, junction densities (P-P, P-R, and R-R 
Junc.D), determined by counting junctions of 1 -voxel-thick 
trabecular structure, did not correlate between the two 
modalities. Therefore, conclusions regarding these parameters 



should be interpreted cautiously. Interestingly, however, Junc.D 
of HR-pQCT significantly contributed to the estimated elastic 
moduli of trabecular bone. Furthermore, they were highly 
discriminative in detecting differences between premenopausal 
women with and without lOP. 

The current HR-pQCT and ITS-based morphologic analysis has 
several limitations. A global threshold technique provided by the 
HR-pQCT manufacturer was used to obtain the input image for 
ITS-based analysis. This threshold technique has the tendency to 
overestimate trabecular bone volume, number, and thickness 
measured by model-independent methods. In addition, the 
image voxel size, 82 (xm, is close to the low end of trabecular 
thickness. A subvoxel process'"*^' followed by a local adaptive 
threshold technique''*''' would be expected to improve the voxel- 
based microstructural and mechanical measurements of images 
in the low-resolution regime. In future studies it will be of interest 
to apply these techniques to in vivo HR-pQCT images to improve 
voxel-based analysis. We expect that these preprocessing 
techniques will also improve the correlations of the ITS 
measurements of HR-pQCT with their gold standard jjiCT 
measurements. 

In conclusion, ITS-based morphological analysis of HR-pQCT 

images detected significant reductions in trabecular plate and 
rod number and bone mass, leading to a less connected, less 
longitudinally aligned, and more widely separated trabecular 
network and reduced bone mass at the distal radius and distal 
tibia in premenopausal women with unexplained osteoporosis. 
On the other hand, trabecular plate and rod thickness remained 
unchanged. A transition to a more rod-like trabecular micro- 
structure was found in the distal radius but not in the distal tibia. 
Most ITS measurements contributed significantly to the elastic 
moduli of trabecular bone over and above the bone volume 
fraction. This study provides in vivo evidence that trabecular 
plates play a far more important role than rods in determining 
the elastic moduli of trabecular bone. At a fixed bone volume 
fraction, plate-like trabeculae contribute positively to the 
mechanical properties of trabecular bone. These findings 
suggest that ITS-based morphological analysis of HR-pQCT 
images provides additional information regarding bone strength 
and is a sensitive and promising clinical tool for the investigation 
of trabecular bone microstructure. 
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Glossary 



Definition of the ITS-Based Microstructural Parameters 
pBV/TV Plate bone volume fraction: the total 

volume of trabecular plates divided by 

the bulk volume 
rBV/TV Rod bone volume fraction: the total 

volume of trabecular rods divided by 

the bulk volume 
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aBV/TV 

pBV/BV 
rBV/BV 

pTb.N (1/mm) 

rTb.N (1/mm) 

pTb.Th (mm) 
rTb.Th (mm) 
pTb.S (mm'^) 

rTbi (mm) 

R-R Junc.D (1/mm^) 

P-R Junc.D (1/mm^) 

P-P Junc.D (1/mm^) 



Axial bone volume fraction: the total 

volume of trabeculae aligned with the 
longitudinal axis divided by the bulk 
volume 

Plate tissue fraction: the total volume 
of trabecular plates divided by the 
total bone volume 

Rod tissue fraction: the total volume of 
trabecular rods divided by the total 
bone volume 

Trabecular plate number density: 
the cubic root of the total number of 
trabecular plates divided by the bulk 
volume 

Trabecular rod number density: the 
cubic root of the total number of 
trabecular rods divided by the bulk 
volume 

Mean trabecular plate thickness: the 

average thickness of trabecular plates 
Mean trabecular rod thickness: the 
average diameter of trabecular rods 
Mean trabecular plate surface area: 
the average surface area of trabecular 
plates 

Mean trabecular rod length: the aver- 
age length of trabecular rods 
Rod-rod junction density: the total 
number of junctions between rod 
and rod divided by the bulk volume 
Plate-rod junction density: the total 
number of junctions between plate 
and rod divided by the bulk volume 
Plate-plate junction density: the total 
number of junctions between plate 
and plate divided by the bulk volume 
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